Topological phases for bound states moving in a finite volume 
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We show that bound states moving in a finite periodic volume have an energy correction which 
is topological in origin and universal in character. The topological volume corrections contain 
information about the number and mass of the constituents of the bound states. These results have 
broad applications to lattice calculations involving nucleons, nuclei, hadronic molecules, and cold 
atoms. We illustrate and verify the analytical results with several numerical lattice calculations. 
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Over two decades ago Liischer derived a relation con- 
necting the energy levels of an interacting two-body 
system in a periodic cube to physical scattering phase 
shifts [l|, y . This finite- volume technique has become a 
standard tool in lattice quantum chromo dynamics [J-Q 
and in lattice effective field theory for nucleons and cold 
atomic systems |9l4l3|. In this letter we consider finite- 
volume effects of composite particles in motion. We dis- 
cuss corrections to the binding energies of bound states 
in a moving frame. We also show how the finite- volume 
scattering method is modified if one or both particles are 
composite. 

We find topological phase corrections associated with 
the motion of bound states in a periodic box. These 
corrections have a universal dependence on momentum 
determined by the number and mass of the constituents. 
In asymptotically large volumes the corrections are ex- 
ponentially small and can be neglected. However, it is 
often the case in large-scale lattice simulations that cal- 
culations are performed at volumes which are not asymp- 
totically large. Fortunately, we find that the corrections 
have a simple form which can be subtracted out from the 
analysis. 

We start with another result derived by Liischer 
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It describes finite-volume corrections to the binding en- 
ergy of two-body dimer states for interactions with finite 
range. For a dimer at rest, the shift in the energy when 
placed in a periodic cube of volume L^ is 

AE^{L) « ^ fd'r 0:,(r) V{r) 0o,(r + nL) . (1) 

Here V{r) is the interaction potential and (j)oo is the 
infinite-volume wavefunction as a function of the rela- 
tive separation r. The summation is over integer vectors 
ft with magnitude 1. Throughout our discussion, we as- 
sume that the energies and momenta are non-relativistic. 
For finite-range interactions Eq. ([T|) gives a correction 
which scales as e~'^^ / L in the large volume limit, where 
K is the binding momentum. 



For general A^-body bound states a straightforward 
generalization of Liischer's result yields 



Ai^ow- E /[n/ 



Ari,---) 



X V(fi,---)(l)oo{fi+niL,---). (2) 

Here f*i are A^ — 1 relative coordinates and fii are again 
integer vectors. Although Eq. ^ clearly does not apply 
to relativistic quarks held by confinement within a single 
meson or baryon, these corrections are useful for analyz- 
ing lattice calculations of hadronic molecules, nuclei, and 
cold atomic bound states 
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We now consider a dimer moving in the same periodic 
cube with momentum 27rfc/L for integer k. In the dimer 
wavefunction we can factorize out the phase dependence 
due to the center-of-mass motion. 



V-L {ri.r2) = e»2^"'=-^"i/-^e*2'^'^-")'^-*'^/^0L(ri 



r2 



(3) 



where a — mi/{mi + 7712). Since ipL (n.,'^2) is periodic 
in fi and f2 , 0l gets a nontrivial phase for each winding 
around the toroidal topology of the periodic cube. 



(l)L{f+ nL) = e 



-i27:ak-n 



H{r)^ 



(4) 



for all integer n. We note that phase factors have been 
previously discussed in connection with finite-volume 
scattering in moving frames 19l421|. However, the phase 
factors have a qualitatively different effect on bound 
state wavefunctions which simultaneously touch all wall 
boundaries. Each phase twist induces a measurable shift 
in the binding energy. 

When Eq. ^ is combined with Eq. ([2]), we find that 
the finite- volume correction is a sum of sinusoidal func- 
tions of momentum. For S'-wave dimers with momentum 
lixhjL, the finite-volume correction has the form 



A£;g(L) - 3 



\Y. 



cos {2n aki) = T{k, a) . (5) 



1=1,2,3 



Finite- volume corrections for higher angular momentum 
bound states at rest have recently been discussed in [22] ■ 
From the results presented there, it is straightforward 
to derive analogous results for T{k,a) for dimers with 
angular momentum. 

For bound states with more than two particles, the 
finite- volume correction has the same general form. For 
an A^-body bound state with all equal masses and no 
cluster substructure, the topological phase is 



i2-nk-ni/N ^^ / 



4>Li--- yfi+fiiL,---) = e 
For S'-wave bound states we get 



L[--- ,r.,, 
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For A^-body bound states with a two-cluster substruc- 
ture, one can apply the two-body result, Eq. ([5]), with 
rrii and 7TI2 being the masses of the two clusters. For 
more complicated A^-body bound states with more than 
two clusters and/or particles with unequal masses, the 
same cosine functions as in Eq. ([5]) are also expected 
for AEj:{L)/AEj^{L). In these cases, however, more in- 
formation is needed regarding which particles or clus- 
ters occupy the tail of the bound state wavefunction. If 
this is unknown, then AEj:{L) / AE^^{L) can be extracted 
from numerical calculations, and an empirical fit to cosine 
functions as in Eq. ([S]) can yield structural information 
about the tail of the bound state. 

Our results presented above have already been adapted 
by Davoudi and Savage into a general method for reduc- 
ing finite- volume errors for two-body bound states such 
as the deuteron using boosted frames [2J|. The compu- 
tational advantage of this approach is that finite- volume 
effects can be directly removed from lattice data without 
extrapolating to large lattice volumes. This is especially 
useful for the case with more than two constituents where 
the analytic form for the finite-volume L-dependence is 
a priori unknown. 

To illustrate the utility of the boosted-frame method, 
we consider lattice calculations of the triton at lead- 
ing order in pionless effective field theory [23[. We use 
the leading-order lattice action defined in Eq. (4.6) of 
Ref. [10] with spatial lattice spacing 1.97 fm and tem- 
poral lattice spacing 1.32 fm/c. The two-body contact 
interactions, C and C/2 , are set to reproduce the phys- 



ical neutron-proton scattering lengths, a^So 



-23.7 fm 



and 035^ — 5.4 fm. The three-body contact interaction, 
-D, is determined by the triton energy at infinite volume, 
-8.48 MeV. 

Using the Lanczos algorithm for sparse-matrix eigen- 
vector iteration 25], we have computed the triton en- 
ergy as a function of periodic cube length L in spatial 
lattice units and momentum 2'Kk/L. In Table U we show 
a comparison of lattice results for AEj:{L) / AE^{L) for 
the triton versus T(fc, 1/3) for lattice sizes L = 6,7, 8. As 
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L = 6 


L = 7 


L = 8 


r(fe,l/3) 


(1,0,0) 


0.395 


0.432 


0.458 


0.500 


(1,1,0) 


-0.035 


-0.025 


-0.016 


0.000 


(1,1,1) 


-0.376 


-0.413 


-0.442 


-0.500 



TABLE I: Comparison of triton lattice results for 
A£'g(L)/A£'g(L) and r(fc, 1/3) versus L. 
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L = 5 


L = 6 


L = 7 


L = 8 


(0,0,0) 
(1,1,0) 


-0.603 
0.029 


-0.169 
0.006 


-0.049 
0.001 


-0.015 
0.0002 



TABLE II: Triton finite-volume energy corrections AEj:{L) 
in MeV for k = (0, 0, 0) and k = (1, 1, 0) versus L. 



seen in Table U the lattice results approach T{k, i/s) in 
the large-L limit. We note that the leading finite- volume 
corrections vanish for k = (1, 1,0). Therefore the calcu- 
lation of the triton binding energy in this boosted frame 
should converge much more quickly in the limit of large 
L. In Table HIl we show the triton finite- volume energy 
corrections AE^{L) in MeV versus L for k = (0,0,0) 
and k = (1,1,0). We see that the finite-volume errors 
are reduced dramatically for fc = (1, 1, 0). 

We now turn our attention to the scattering of com- 
posite states in a finite periodic cube. We consider the 
scattering between states A and B in the center-of-mass 
frame. Let fiAB be the reduced mass, and let Eab{p, L) 
be the total energy of the A — B scattering system with 
radial momentum p in a periodic cube of length L. States 
A and B can be point particles or composite bound 
states. We assume that the constituent particles com- 
prising the states have finite range interactions. The 
composite structures of A and B, however, will in gen- 
eral produce effective interactions with exponential tails 
extending to infinity. 

These tails generate exponentially-small finite- volume 
corrections to Eab{Pi L) associated with the binding en- 
ergies of A and B separately as well as the scattering of A 
and B together. In this analysis, we focus only on the ex- 
ponential corrections to the binding energies. This will be 
useful in lattice simulations where the volume is not very 
large and the binding energy shifts may be comparable 
to that of the scattering process being measured. We will 
not be concerned with exponentially small corrections to 
the scattering of A and B. If the interactions between 
A and B are very strong, then it is theoretically possible 
that the finite-volume scattering corrections we neglect 
are comparable to the binding energy shifts. However, 
in such cases the part of the energy shift due to scatter- 
ing which is not exponentially suppressed will be much 
larger still, and so the loss of accuracy in the scattering 
analysis will be small. 

In order to calculate finite-volume corrections due to 



the binding energy, it suffices to consider singular solu- 
tions of the free Helmholtz equation. Let f be the sep- 
aration between the center of masses of the two states. 
In the following we assume that p is sufficiently small so 
that angular momentum mixing with higher-order sin- 
gular solutions can be neglected. For S'-wave scattering 
between states A and B with radial momentum p, the 
position-space scattering wavefunction is 



\'^p)-cj2- 



'2Trk/Ly 



(8) 



with some normalization constant c. 

We let E^{L) and E^JL) be the finite- volume ener- 
gies due to binding for bound states A and B with mo- 
menta 27rfc/i and — 27rfc/L, respectively. For point parti- 
cles without internal structure, these energies are by def- 
inition zero for all momenta. The total energy Eab{p^ L) 
is then 



Eab{p,L) 
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(9) 

where N = Y.k (^^ - v) ^ and t? = p'^ L"^ j {luf . The 
finite- volume correction can be written as 

Eab{p, L) - Eab{p, oo) 

= TAiv) AE^{L) + TBiv) ^E?{L) , (10) 

where AE^{L) and AE? {L) are the finite- volume cor- 
rections for states A and B at rest, and we have defined 
the topological volume factor 

1 Y^ 121 = 1,2,3 '^°s(27rafci) 



t{7j) = 



k 



3{P-nj 



(11) 



The analysis can be generalized to higher angular mo- 
mentum scattering states using an extension of Liischer's 
scattering relation to higher orbital angular momen- 



26|,l27|. 



tum 

The finite-volume correction in Eq. (jlOp has nothing 
to do with the interaction between states A and B and 
should therefore be subtracted from the total energy be- 
fore using Liischer's scattering relation. This subtraction 
should reduce systematic errors in lattice calculations. 

As an example to test the method, we consider 
fermion-dimer scattering for two-component fermions. 
The physical process which we study corresponds with 
spin-quartet scattering between a neutron and deuteron. 
In that case the two fermion components should be re- 
garded as protons and neutrons with the same spin 

We use the same lattice Hamiltonians as in Ref. 
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except in that case the scattering length was tuned to 
infinity. As in Ref. [28|, we consider two different lattice 
Hamiltonians, each of which produces a shallow dimer 



in the continuum limit. The first Hamiltonian, Hi, con- 
tains only a local contact interaction between fermions. 
The second Hamiltonian, H2 , contains a contact interac- 
tion as well as nearest-neighbor interactions. These are 
used to tune the binding energy of the dimer while also 
setting the effective range parameter to zero. Both lat- 
tice Hamiltonians reproduce the same continuum limit of 
fermions with attractive zero-range interactions. 

We now focus on the fermion-dimer system. This cor- 
responds with a neutron together with a deuteron in the 
spin-quartet channel in pionless effective field theory at 
leading order. Experimental measurements find a quar- 
tet scattering length '^Und = 6.35(2) fm [29|. This can 
be expressed as a fraction of the spin-triplet neutron- 
proton scattering length, ^a^dl'^anp = 1.17(1). A more 
detailed calculation includin g in teraction range effects 
obtains ^a„d = 6.33(10) fm [30|, |3l[, in full agreement 



with experimental values. 

We calculate fermion-dimer scattering on the lattice 
using Liischer's finite-volume formula. Using the Lanc- 
zos algorithm we have computed the ground state energy 
for the fermion-dimer system on periodic cubes for six 
different lattice spacings aiatt- For each lattice spacing 
we consider periodic volumes I? ranging from L = 6 to 
L — YJ lattice units. From these energies we determine 
the low-energy parameters for fermion-dimer scattering 
and extrapolate to the continuum limit. The full details 
of this calculation will be described in a forthcoming pub- 
lication, and we just summarize the results here. 

Results for the fermion-dimer scattering length, apD, 
and the effective range parameter, rpD, are shown in 
Fig. [1] In each case, we extrapolate to the continuum 
limit and write final results as dimensionless combina- 
tions multiplied by powers of the dimer binding momen- 
tum, kd. In the shallow binding limit kd equals the recip- 
rocal of the fermion- fermion scattering length. To see the 
effect of the topological volume correction we have done 
the full calculation using the correct topological factor 
T(r\), as well as a faulty calculation which replaces T[r\) 
by 1. For comparison, we show the continuum results, 
apD^D = 1.17907(1) and tfdKd = —0.0383(3), obtained 
via the Skorniakov-Ter-Martirosian (STM) integral equa- 



tion 
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33| . We find that the small size of the effective 
range parameter requires a fit to low-energy scattering 
that includes the shape parameter, which was not done in 
earlier calculations of the effective range parameter 34 1 . 
We note also the agreement with the experimental value 
^a„d/'a„p = 1.17(1). 

The results in Fig. [1] show that the inclusion of the 
topological volume factor t(j]) is essential for obtaining 
the correct continuum limit. In all cases the continuum 
extrapolations for B\ and Hi agree within error bars. 
However, the correct T(r\) factor is needed to reproduce 
the STM equation result. The correction is small for the 
scattering length, but quite large for the effective range 
parameter. 
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FIG. 1: (Color online) Lattice results and continuum extrap- 
olation with error estimates for the fermion-dimer scatter- 
ing length (top) and effective range parameter (bottom). For 
comparison we show the continuum results obtained via the 
Skorniakov-Ter-Martirosian equation. 

We expect the topological volume factor to have im- 
portant effects in other lattice calculations of scattering 
for composite bound states. The analysis presented here 
should provide a simple but general method for improv- 
ing the accuracy of bound-state scattering calculations. 
The list of possible applications is quite extensive and 
include lattice calculations involving the scattering of 
nucleons upon nuclei, the scattering of nuclei, Compton 
scattering and electroweak probes upon nuclei, mesonic 
and baryonic scattering upon hadronic molecules, and 
few-body scattering in cold atomic systems. An appli- 
cation of this method to dimer-dimer scattering is in 
progress. 
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